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Abstract 

It is shown that a strongly magnetized isothermal pair plasma near the surface of a pulsar supports 
low-frequency (in comparison to electron cyclotron frequency) toroidal electrostatic plasma modes 
in the equatorial region. Physically, the thermal pressure coupled with the magnetic pressure 
creates the low frequency oscillations which may grow for particular case of inhomogeneities of the 
equilibrium magnetic field and the pair plasma density. 
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1 Introduction 

Study of pair plasma in the pulsar magnetosphere 
is related to the investigation of radio emission coming 
from these sources (Buzzi and Hines, 1995; Mofiz, 1997 
and the references therein). Radio pulsars, which are 
rotating neutron stars with spin periods ranging from 
1.57 ms to 5 ms, are characterized by surface magnetic 
fields in the order of lO^^G, radii of about 10 km. 
A rotating magnetized neutron star generates huge 
potential difference between different parts of its surface 
(Goldreich and Julian, 1969). Theoretical models (Stur- 
rock and Baker 1979; Ruderman and Sutherland, 1975; 
Arons and Scharlcmann, 1979) have been developed to 
predict the production of pair plasma (electron-positron) 
in the pulsar magnetosphere. A mechanism of cascade 
generation of particles has been suggested, according to 
which secondary electrons and positrons result from pair 
production induced by high-energy curvature radiation 
photons- emitted by primary electron beams coming 
from pulsar surface. Pair plasmas are composed of 
charged particles with same mass and opposite charges, 
which admit the time and space scales that significantly 
differs from those of an electron - ion plasma . There 
appear a great variety of linear and nonlinear plasma 
modes in the strongly magnetized pair plasma. A good 
number of works are done by several authors ( Sakai and 
Kawata, 1980; Lominadze et al., 1983; Mofiz, 1989, 1990, 
1997;Mofiz et aZ.1985, 1987; Shukla et al.,2007). Recent 
discoveries (Sana et a/. ,2010) of kilohcrtz quasi-periodic 
oscillations (kHz QPOs) may also be related with the 
typical dynamical time scales of pair plasma close to the 
surface of the neutron star. 



In this paper, we have investigated the low frequency 
(in comparison to the electron cyclotron frequency) 
electrostatic oscillations in the equatorial region of the 
pulsar magnetosphere. We found that thermal pressure 
coupled with the magnetic pressure creates the low fre- 
quency oscillations which propagate along the toroidal 
direction of the equatorial plasma. The magnetic field 
inhomogeneities and the equilibrium plasma density in 
particular situations may admit growing low frequency 
electrostatic modes in the pulsar magnetosphere. 

2 Formalism: Equilibrium Magnetosphere 

A. Equilibrium magnetic field 

To have an idea of the equilibrium plasma density in the 
pulsar magnetosphere, we consider the simple model of 
dipolar nonalighned magnetic field B — {Br, Bg,0} of 
the pulsar (Ginzburg and Ozernoy, 1964; Muslimov and 
Harding, 1997; Mofiz and Ahmedov,2000; Prasanna et. 
al., 1989 ) with the components: 
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where, Br, 



is the Newtonian value of the magnetic 
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field at the pole of the star with the magnetic moment 
fi and radius R; rg is the Schwarzschild radius of 
the pulsar( Shapiro and Teukolsky,2007).The detail 
calculation of magnetic field structure around a 
nonaligned neutron star was given by Ginzburg 
and Ozernoy in 1964, which was further devel- 
oped by different authors including Muslimov 
and Harding in 1997. Here, we consider the 
expression for the components of magnetic field 
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For Newtonian case (no gravity) e = 0, we have 

f = Csin^e, (7) 

Now, in order to visualize the field line structure, it 
is useful to transform over to a Cartesian frame through 
the usual relation ( X = rsinOcos^, Y — rsinOsin^, Z = 
rcosO). Then, from Eq.(7) the magnetic lines of force 
projected on the meridional plane ($ = 0) are described 
by the equation 
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FIG. 1: Magnetic field at the pole and at the equatorial region 
of the pulsar magnetosphere : solid line ( — ) at the pole; 
broken line (— —) and doted line (...) at the equator for 
e = 0, .5, respectively. 



Br,Bg, which are approximately valid near the 
surface of a pulsar 



Introducing a compactness parameter e — 

Bojrfi) 



R ■ 



and dimen- 
we can write the magni- 
tude of the pulsar magnetic field as: 



mensinless magnetic field Bo{r,9) — g 
sionless radial distance r 



-Bo(r,( 



co8^e+ -(I - ^)sin'^9. (3) 
4 r 



The magnitudes of the magnetic field at the pole and 
at the equatorial region of the pulsar magnetosphere, 
are shown in Fig.[l] The figure shows that the magnetic 
field is more strong at the pole , where gravity has no 
effect and it is less at the equator where gravity reduces 
the field. 

To see the magnetic field configuration, we solve the 
equation for magnetic lines of force (Banarjee et al, 1995) 



dr rdO 
Br Ba ' 



(4) 



which gives 
2r 



(^l + ^1-^^ -rg^ 4Cexp2(V^-i) s^n^0^ (5) 

where C is a constant of integration. 

Eq.(5) can be written in dimensionless form as 



where, x = = j^. 



For compact object with gravity (e ^ 0), the equation 
for lines of force [ Eq.(6)] can be written as 



2(1 + yr^) 

which for 4^1, can be written as 



(9) 



_3 ^ 4Cx^exp^(V^-^)+efg ^^^^ 

2(1 +yr^) 

where, we put ^ ~ ^r- with = \J Cx^ . 

Then, the magnetic lines of force in the meridional 
plane are described by the equation 



4Cx2 exp^^V^ 



2/3 



(11) 



Using Eqs. (8) and (10), we plot the lines of forces 
for Newtonian case (e = 0) and for gravitational case 
(e = .5) putting the values of (7 = 1,2,3,5; respectively, 
which are shown in Fig. [2] 



To see the visual observation of the ambient magnetic 
field around a pulsar, we perform the density plot of Bq 
from the Eq.(3) in the meridional plane {x,z) which is 
described by the equation 
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where, we put C 
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The distribution of magnetic field magnitudes are 
shown in Fig. [3] and Fig. [4] for the Newtonian and 
gravitational cases, respectively. 
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FIG. 2: Magnetic lines of force in tiie meridional plane (a;, z) 
of the pulsar magnetosphere:solid line ( — ) [Newtonian case 
e = 0], broken line(— —) [gravitational case,e = .5]. 
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B. Equilibrium plasma density. 

Now, we study the equiUbrium plasma density around 
pulsar. The equilibrium is due to the balance of magnetic 
pressure and the thermal pressure, which is described by 
the equation (Krall and Trivelpiece,1973) 



Stt 



In spherical coordinate system: 
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Then, Eq.(13) yields: 

2no(r,0)To = / ^(^= 

47r J r /i _ Is 
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l)dr 



Now, we investigate the equilibrium equation (16) for 
the following two cases: 



i) Newtonian case 



FIG. 3; Distribution of magnetic field magnitudes 
in the meridional plane (a;, 2)of the pulsar magneto- 
sphere: (Newtonian case, e = 0). 
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FIG. 4; Distribution of magnetic field magnitudes 
in the meridional plane (a;, z)oi the pulsar magneto- 
sphere: (gravitational casee = .5). 



For a Newtonian case, = 0, so 
B,(r,0) =Bp(|)'cos0, Bg{r,e) 
Then from Eq. (16), we have 



Bp {R\^ 
2 



■p) sinO 



where Ci(9) and C2(r) are the constants of integra- 
tions, which are to be determined from the boundary 
conditions. 



We set the boundary conditions: i)at r oo; i?o ^ 
and no — > 0, ii) at r = i? and = 0; ng — > rip, where. 
Up is the equilibrium plasma density at the pole of the 
pulsar. Then, we find Ci{9) — and C2(r) = const — 
-47r(l + Thus, from Eq.(17), we get 



R 

no{r,9) = no{r) = np{ — f, 



(18) 



Equlibrium plasma density (Newtonian case) 




FIG. 5: Equilibrium plasma density no (r) (Newtonian case, 
. = 0) 
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FIG. 6: Equilibrium plasma density distribution in the merid- 
ional plane {x, z) (Newtonian case,e = ) 



Then, from Eq.(16) we find 



877 



2no(r,0)To 



Idtt r 

4 .6 



Li. 
7r 



9009 V, 



"a 

Bl,R 



r r 



47r r 



r 



1 - 



+ D2{r) 



4r(l-^) 

(20) 



where Di{9) and D2{r) are the constants of integration 
and 
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The constants Di{9) and D2{r) may be determined 
from the earlier boundary conditions, i.e.:i) at r — > oo; 
Bq ^ and no — and ii) at r = i? and 6 = 0; hq = rip, 
where is the equilibrimn plasma density at the pole. 
Thus, we find 



Di{9) = const. 
£'2(7') = const. = - 



4 X 256 
9009 
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(22 
2rtpTo 



Then, from Eq.(20), we obtain the equilibrium plasma 
density under gravity as 



which is independent of 6 and can be represented as 

,1 



no(r-) = {-_f 
r 



(19) 



where, no = It shows that in the Newtonian case, 
the plasma density around the pulsar is isotropic and is 
a steeply decreasing function of radial distance, which is 
shown in Fig. [5]. The corresponding density distribution 
in the meridional plane {x,z) is shown in Fig. [6]. 



ii) Gravitational case 



In the case of gravity ^ 0, so 
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Equlibrium plasma density at the equator (gravitational case) 
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FIG. 7: Equilibrium plasma density at the equatorial disc 
of the pulsar magnetosphere under gravity: solid line ( — )for 
e — .7; broken line ( )for e — .5; dot line (...)for e = .1 



which can be written in dimensionless form as 



ujc = eBo/meC, where e is the charge of electron, me 
is the mass of electron, c is the speed of light ) elec- 
trostatic oscillations in the equatorial region {9 = n/2). 
In this case, the perpendicular ( to ^ ) components 
of electron/positron velocities in the electrostatic field 
E± = — Vj_<I>, where <I> is the electrostatic potential, are: 



Boir) 



X V I $ 



qjBoir)noir) 



9 X V±nji, (26) 



where qj is the charge of species j = (e,p) and nji{<< 
no(r)) are small electron/positron number density per- 
turbations. 

The number density perturbations of the plasma fluids 
are determined by the continuity equation 
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(27) 



which using the Eq. (26) can be written as 
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Here, we introduce the dimensionless quantities: 
no(f,^) = ^and5,^ = 5^. 

Now, for different values of compactness parameter 
e — .1, .5, .7, respectively; we plot the equilibrium plasma 
density in the equatorial disc of the pulsar magnetosphere 
which are shown in Fig. [7] . 

The above figure (Fig. [7]) shows that the equilib- 
rium plasma density varies significantly for increasing 
values of compactness parameter e. Very close to the 
surface of the compact object, for a higher value of the 
compactness parameter e = .7, the equilibrium plasma 
density sharply increases but away from the surface it 
falls steeply as in the Newtonian case. 



3 Electrostatic modes in the equatorial region 



Against the equilibrium, now we study the low fre- 
quency (in comparison to electron cyclotron frequency 



driji 



1 d 



cno{r) 



dr \ Bo{r) 
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Since, in the present study we are interested only in the 
toroidal modes , we consider rt^iand $ to be proportional 
to exp (—iwt -|- jmi/)), where m is an integer 0, ±1, ±2,.... 
which represents the toroidal mode number and u is the 
frequency of the electrostatic mode. Then, from Eq.(28), 
we get 



m „_§_(„ rio(r)_\ 
_ r^^dry Bo(r)l 
"j'l ~ m d /- cT, 



^ 7^ drV q,Bo{r)) 

From the Poisson's equation 

Vi$ = 47re(nei-npi), 
we find the dispersion relation 
^2 /I a / cT 



(29) 



(30) 



UJ 



j,2 y-r dr \ eB^ir) 
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where, we put — Tp ~ T ~ To/2. 



Now defining 

19 A ^ i^fr^) = Vb with Vt = 

r or \ eBolr) J r or \ uJc J ^ ^ 

/\1/'^ 1/2 

I J and defining ujp — (Sne'^no/me) for electron 
positron plasma, we can rewrite the dispersion relation 
(31) as 

^2 

uj^ - KM + -^KnbVB = 0, (32) 
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where Kl = ^ and K„6 = -^-4-{r'^). 
Letting a; = cjr + ir in Eq.(32), we get 

^2 

w2 _ r2 - KlVl + -^Kr^^VB - 0, 

2tJrr = 0, 



(33) 



(34) 



Eq.(34) gives two cases, namely, either a)^;^ ~ 
or b)r — 0. For cj,. = from Eq. (33), wc obtain the 
growth rate 
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Combined effect of inhomogeneides for e=.7 




which shows the electrostatic modes grow if 



i^nh > 2~- 

UJr, 



(35) 



(36) 



For a detail investigation, let us consider the o-mode i.e. 
m — 0. In this case — 0, then, we find 



c / no{r)T 
ry BUr)/87r 



where 



A(r) = (— - — -1)(1-— ), 

J^B J^n J^B 



(37) 



(38) 



with Lb = ( 



1 O-Bq- 
Bn 



and Lr. 
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1 dno^ 



representing 



^ ^ -tin Or ' ^ no c)r ^ ^ 

the scale lengths of inhomogeneities of the magnetic field 
and the equilibrium plasma density, respectively. In fact 
, A(r) represents the combined effect of inhomogeneities. 



In Eq.(37), represents the ratio of thermal pres- 

sure to the magnetic pressure, which is a positive quan- 
tity depending on r. Therefore, for F > 0, A(r) is to 
be positive. The graphical representation of X(r) for a 
particular value of e = .7 is shown in Fig. [8] which shows 
that A(r) is positive for approximately r > l.li? , it is 
negative for r < I.IR and for r ^ i?, it is saturated. 
The saturation is due to the fact that for r ^ R, 
j-^ ^ —i-ij^ ^ —6 and as a whole A(r) ~ 8 which is 
a constant. Thus in the case of high compactness of a 
pulsar, the generated electrostatic oscillations in the pair 
plasma at the equatorial region may grow. For F = 0, 
from Eq. (33), we obtain the frequency of oscil- 
lations which is given by 



c / no(r)T 



(1 



r 



-)(1 



-)• (39) 



FIG. 8: Combined effects inhomogeneities of magnetic field 
and plasma equilibrium density at the equatorial disc of the 
pulsar magnetosphere for the compactness parameter e = 0.7. 



which shows that u 



c( no(r)T n1/2 
r yBl(r)/?,TTl 



, where no{r)T 

is the thermal pressure and B'^{r)/d>Ti is the mag- 
netic pressure. Thus, it seems that the coupling 
between the thermal pressure and the magnetic 
pressure, which creates the electrostatic oscilla- 
tions in the pulsar magnetosphere. Although the 
frequency of oscillations is much lower in com- 
parison to electron cyclotron frequency which is 
extremely high for a pulsar magnetosphere, in 
fact, a reasonably high frequency electrostatic os- 
cillations may be created through this coupling. 
This coupling also may create high frequency 
electromagnetic oscillations, which needs further 
scrutiny. 



At r 



R, we find lu 



no(R)T 



For 



R\l B^i?)/87r- 

a pulsar with super strong magnetic field, we may 
consider that magnetic pressure is much higher than 
the thermal pressure i.e. Bq{R)/8tt ^ no(i?)r near the 
surface, then ujr For crab pulsar PSR 0531 +21, 

Bo{R) = lO^^G and R = lO^cm (Nanobashvih, 2004). 
So we find w_r < 3 x lO^^s^^ and ~ lO^Os^^, which 
validate our approximation w <C and accordingly 
we ignored the effects of currents and conduc- 
tion. 

So far, we consider the nonaligned magnetic 
field. It would be interesting to investigate 
the growing modes in oblique direction of the 
magnetic field. To find clues in this case, we have 
investigated the simple case of aligned magnetic 
field (Michel and Li, 1999), which shows the 
mode grows from a particular range of angle 

inclination a > sin^^ ^ ^ _b" /"stt • 

4 Discussion 

Here, we study the low frequency (in comparison to 
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the cyclotron frequency)electrostatic modes in the 
isothermal electron-positron plasma of the pulsar mag- 
netosphere. The compactness of the pulsar is considered 
in this investigation. The equilibrium magnetic field and 
the equilibrium plasma density are investigated in detail. 
To visualize the picture of the magnetic field around 
pulsar, the magnetic lines of force are drawn for both the 
Newtonian case (no gravity) and for the compactness 
case (with gravity). The corresponding magnetic field 
distribution in the meridional plane are shown in density 
plots. The equilibrium is due to the balance of thermal 
pressure with the magnetic pressure. In the Newtonian 
case, it is found that the equilibrium plasma density 
is isotropic and is steeply decreasing function of radial 
distance. But for the compact pulsar with gravity, 
the equilibrium plasma density varies significantly for 
increasing values of compactness parameter. Further, 
against the equilibrium, we study plasma perturbation 
in the linear approximation. The number density 
perturbations of the plasma fluids arc determined 
by the continuity equation, which coupled with the 
Poisson's equation, yields the dispersion relation for the 
electrostatic modes. From the analysis of the dispersion 
relation, we find growing modes, which are related 
with the inhomogcneity scales of equilibrium magnetic 
field and the equilibrium plasma density. The thermal 
pressure coupled with the magnetic pressure creates the 
low frequency electrostatic modes which propagate along 
the toroidal direction of the pulsar. 

5 Conclusion 

Growing elecrostatic modes may be generated in the 
equatorial region of a compact pulsar magnetosphere. 
The modes are created due to the coupling of the 
thermal pressure with the magnetic pressure, which may 
grow for the particular cases of inhomogeneities of the 
equilibrium magnetic field and the equilibrium plasma 
density. This may have some relevance to the recent 
discoveries of low frequency quasi-periodic oscillations 
(QPOs) of magneter with an extremely strong magnetic 
field (Lee, 2008). 
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